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LABORATORY  PROCEDURES 


ION  AUGMENTATION  EXPERIMENTS 


l^he  conduct  of  the  experimental  analysis  of  electron  seeding  techniaa 
uig  lUMuiiLdUy  Ouii'lmuHKu.  iW'  qwpwpglg^ill^^ 

to  clarify  the  Experimental  technique  s^T^cT^to  detail  th^  mathematics^ 


procedures  necessary  for  analysis  and  interpretation  of  the  data# 

The  basic  purpose  of  the  experiments  is  twofold:  first,  quantitative 
empirical  definition  of  the  performance  of  electron  augmenters,  both  ablative 
and  rocket  motor  types,  when  subjected  to  a  range  of  conditions  which  will 
affect  their  performance;  second,  an  estimate  of  the  important  processes 
and  their  rates  which  determine  the  rate  of  electron  disappearance  in  flow 
subsequent  to  their  generation.^ 

2.  Planned  Experiments  ^  p  oJ' ^ 


Rocket  motor  augmentation,  ablation  augmentation,  and  chemical 
kinetics  experiments  are  currently  planned.  The  first  two  sets  of  experiments 
are  called  for  in  the  ciarrent  Bendix  Program  Plan.  The  third,  while  not 
specifically  required,  is  considered  necessary  for  an  accurate  estimate  of 
electron  depletion  rates  due  to  a  paucity  of  available  data. 


2.  1  Rocket  Motor  Experiment 

The  basic  purpose  of  the  rocket  motor  experiment  is  to  determine 
the  empirical  electron  flow  rate  from  the  nozzle  exit  plane  for  a  specific 
rocket  motor  design. 


2.  1.  1  Vacuum  Chamber  Operation 

The  rocket  motor  will  be  operated  in  the  vacuum  system  as  shown 
in  Figure  1  at  exhaust  pressures  corresponding  to  altitudes  of  100,  175, 
and  225  kilofeet  altitude  (8.  29  mm  Hg.  ,  0.  443  mm  Hg.  ,  0.  0586  mm  Hg.  , 
respectively).  These  pressures  are  nominal  settings  prior  to  the  run. 
Variations  of  -  10%  are  acceptable.  The  actual  pressure  during  rocket 
operation  is  to  be  recorded  continuously. 
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Figure  1:  General  Equipment  Layout 


2.  1.2  Rocket  Motor  Ballistics 

The  rocket  motor  will  be  lu  n  according  to  the  procedures  set  up 
by  the  Atlantic  Research  Corp.  Pressure  and  temperature  will  be  recorded 
continuously  via  the  pressure  sensors  and  thermocouples  provided.  The 
weights  of  active  components  (propellant,  igniter,  and  insulator)  shall  be 
measured  before  and  after  each  run  and  the  average  mass  flow  of  each 
determined  from  the  running  time. 

For  each  altitude  the  nozzle  shall  be  nominally  correctly  expanded 
based  on  estimated  frozen  gas  properties  and  nominal  chamber  pressure. 
Although  it  would  be  desirable  to  have  uniform  parallel  flow  at  the  nozzle 
exit,  the  uncertainties  in  the  pressure  and  various  gas  properties  do  not 
make  it  practical  to  design  such  a  nozzle.  The  specific  nozzle  geometries 
are  shown  in  Figure  2. 

2.  1.3  Spectrographic  Measurements 

All  of  the  thermodynamic  and  physical  measurements  made  on  the 
gas  after  it  leaves  the  nozzle  exit  are  spectrographic.  The  layout  of  the 
optical  train  is  shown  in  Figure  1, 

The  optical  system  as  shown  images  the  entrance  slit  of  the  spectro¬ 
graph  to  the  vertical  plane  passing  through  the  center  line  of  the  jet.  The 
vertical  position  of  the  slit  may  be  varied  to  obtain  radial  variation  of 
intensity. 

Since  the  optical  system  is  stigmatic  there  is  a  defined  "resolution 
cell"  from  which  all  radiation  reaching  the  detector  comes.  The  geometry 
of  the  cell  is  analyzed  in  Appendix  A.  The  optical  system  and  the  detector 
of  the  spectrograph  are  calibrated  on  an  absolute  intensity  basis  against  a 
National  Bureau  of  Standards  ribbon  filament  lamp.  The  calibration  pro¬ 
cedures  and  data  are  described  in  Appendix  B. 

The  basic  output  of  the  instrumentation  is  the  radiant  energy  from 
the  resolution  cell  in  ergs  per  volume  of  the  cell  per  spherical  solid  angle. 
The  solid  angle  is  defined  by  the  effective  "f"  number  of  the  optics  per 
wavelength  interval  and  the  effective  wavelength  interval  is  in  turn  deter¬ 
mined  by  the  spectrograph  dispersion  and  effective  slit  width. 
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The  corrected  data  is  of  the  form: 


ergs 


corr 


cm^,  steradian,  X,  sec 


over  the  range  of  wavelengths  and  spatial  extent  desired.  The  fundamental 
equation  for  specie  density  used  in  the  quantitative  spectral  analysis  is 


AX 


I  d  X 
corr 


Ni  gi  A  ij  N  I'.. 


(^) 

sec 


(1) 


where  the  integral  is  taken  over  the  whole  line  of  a  particular  transition 
is  the  number  density  in  the  upper  energy  state  is  the  degeneracy 

of  the  upper  state,  A^j  is  the  transition  probability  from  state  i  to  state  j, 
and  hl^ij  is  the  energy  of  a  single  transition  from  state  i  to  state  j.  The 
number  of  interest  is  Ni  which  is  the  only  unknown  quantity  in  the  equation. 

The  other  factors  are  either  measured  or  known  beforehand.  The  basic 
assumption  made  in  the  derivation  of  Equation  (1)  is  that  self  absorbtion 
is  negligible,  an  assumption  that  improves  as  the  density  is  decreased 
and  is  valid  for  practical  purposes  over  the  range  of  density  to  be  encountered 
in  this  experiment.  The  major  difficulty  expected  in  the  analysis  of  the  data 
is  to  resolve  the  line  of  interest  from  the  background.  In  principle,  the  total 
density  of  a  given  specie  In  a  gas  sample  is  obtained  simply  by  summing  the 
Ni  determined  by  Equation  (1).  In  practice,  however,  not  all  transitions 
are  accessible  to  the  spectrograph  due  either  to  insufficient  intensity  or 
wavelengths  out  of  the  band  pass  of  the  instrument.  If  the  interval  states 
of  a  given  specie  are  in  equilibrium  then  Boltzmann's  distribution  law 
applies  and  the  number  density  may  be  parameterized  in  terms  of  a 
"temperature**  and  the  missing  elements  in  the  sum  determined.  The 
applicable  form  of  Boltzmann's  law  is: 


- 

kT 


N 


(2) 
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The  denominator  of  Equation  (2)  is  the  partition  function,  and  for  the  range 
of  temperatures  considered  here  is  roughly  constant.  In  fact,  little  accuracy 
is  lost  in  calling  it  constant  for  this  calculation.  If  greater  accuracy  is 
necessary,  one  or  two  iterations  will  yield  the  correct  value.  Thus  calling: 


we  may  plot 


i/kT 

=  Q  =  constant  (from  the  experimental 
data) 


Ni  Q  1 

y.  -  -  =  -  T 

&1 


(3) 


versus  G  ^  which  will  be  a  straight  line  with  slope  -  if  the  assumption 
of  internal  equilibrium  is  correct.  The  parameter  T  so  determined  may 
then  be  used  to  calctilate  the  inaccessible  and  the  total  specie  density 

determined  as 


00 

Ns=  X  N3. 

o 

where  '*s”  refers  to  s^^  specie.  In  practice  because  of  the  exponential 
form  of  the  population  equation  the  infinite  sum  is  not  required  (or  possible) 
and  the  sum  may  be  truncated  at  a  few  terms. 

It  should  be  noted  that  the  parameter  T  should  in  no  case  be  inter¬ 
preted  as  the  local  kinetic  temperature  unless  it  is  determined  by  other 
means  that  local  thermodynamic  equilibrium  prevails. 

The  method  described  above  for  specie  density  measurement  is 
quite  general  although  most  conveniently  applied  to  atomic  species  with  well 
separated  lines.  For  molecular  species  it  is  more  convenient  to  deal  with 
integrated  vibrational  band  intensities  of  the  individual  electronic  states. 
Since  Cs  and  Cs*^  are  the  most  important  species  from  the  standpoint  of 
electron  generation  and  depletion,  the  experimental  measurements  will 
be  confined  to  them  unless  experience  shows  that  other  molecular  and 
atomic  species  are  playing  an  important  role. 


6 


The  most  interesting  specie  from  the  standpoint  of  program  ob¬ 
jective  is,  of  course,  the  electrons.  The  preceding  methods  are  not  applic¬ 
able  to  measurement  of  electron  densities  and  other  techniques  must  be 
found. 

The  general  method  involves  the  measurement  of  line  shapes  which 
can,  by  various  theoretical  methods,  be  related  to  the  ion  and/or  electron 
density,  Note  that  it  is  always  assumed  here  that  the  ion  and  electron 
densities  are  equal.  The  details  of  the  analytical  procedures  are  given  in 
Appendices  (C)  and  (D).  Since  the  spectrograph  is  assumed  to  have  been 
calibrated  absolutely,  a  single  spectral  measurement  is  sufficient  for  both 
types  of  analyses,  although  for  the  electron  density  measurements  only  a 
relative  calibration  is  required  and  then  only  to  ensure  that  no  detector 
response  anomalies  exist  in  the  vicinity  of  the  line  in  question.  This  is  so 
since  only  the  relative  line  shape  is  important. 

The  measurements  of  most  significance  are  those  at  the  exit  plane 
of  the  nozzle  since  it  is  only  there  that  there  is  a  reasonably  good  knowledge 
of  the  flow  velocity  and  temperature.  Downstream  of  the  nozzle,  diffusion 
and  mixing  gradually  destroy  the  homogeneous  character  of  the  stream  and 
although  the  measurements  can  be  made  their  interpretation  in  terms  of 
rate  processes  becomes  less  and  less  meaningful. 

A  tentative  selection  of  axial  stations  in  terms  of  nozzle  exit  dia¬ 
meters  (which  will  vary  with  design  exit  pressure)  is  0.0,  0.5,  1.0,  3.0, 
and  5.  0  diameters. 

It  should  be  noted  that  the  0.  0  station  is  nominal.  The  actual 
position  of  the  slit  image  should  be  such  that  no  part  of  the  nozzle  wall  is 
seen  by  the  optical  system  as  shown  in  Figure  3.  This  position  must  be 
found  by  trial  and  error.  The  temporal  sequence  of  the  measurements 
should  be  from  station  0.  0  to  downstream  since  the  limited  running  time 
may  preclude  the  last  measurement.  The  initial  runs  are  to  be  made  with 
the  optical  axis  of  the  instrument  intersecting  the  center  line  of  the  jet 
(Figure  4).  The  analysis  is  to  be  carried  out  assuming  no  radial  gradients 
in  gas  properties.  If  time  permits,  this  assumption  will  be  checked  by 
making  measurements  with  the  slit  displaced  toward  the  edge  of  the  jet. 
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Figure  3  Adjustment  of  Station  0.  0  To  Eliminate 
Nozzle  From  Field  of  View 
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Configuration  for  Initial  Tests 


Displacement 
optical  axis 


Figure  4;  Configurations  for  Tests 


Z.  2  Ablation  Experiments 


The  basic  measurements  to  be  performed  in  the  ablation  testa  are 
similar  in  nature  to  those  of  the  rocket  motor  tests  and  will  not  be  re¬ 
described  here.  The  differences  lie  largely  in  the  experimental  geometry 
and  the  source  of  the  electrons.  The  setup  is  shown  in  Figure  5.  The 
test  gas  (  which  may  be  anything  required  for  the  purpose  of  the  experi¬ 
ment)  is  heated  in  the  arc  chamber,  allowed  to  equilibrate  in  the  setting 
chamber,  and  then  correctly  expanded  to  the  desired  free  stream  conditions 
which  are  the  same  as  in  the  previous  experiment. 

The  ablation  sample,  which  ordinarily  will  be  a  pointed  rod  of  car¬ 
bon  infiltrated  with  a  seeding  material,  is  to  be  placed  a  short  convenient 
distance  from  the  nozzle  exit  plane. 

An  attached  weak  shock  will  develop  at  the  tip  of  the  rod  and  at 
the  shoulder  the  expansion  fan  will  re -expand  the  gas  to  free  stream  pres¬ 
sure.  It  is  assumed  that  because  of  the  weak  shock  the  flow  is  nearly 
isentropic  and  the  free  stream  temperature  also  returns  to  its  initial 
value  (note  that  the  essentially  non  isentropic  flow  can  be  calculated  if 
necessary,  particularly  if  blunt  samples  are  employed).  The  boundary 
layer  will  form  on  the  surface  of  the  rod  and  grow  in  thickness  down  the 
rod.  Conditions  in  the  boundary  layer  are  the  prime  object  of  this  series 
of  measurements. 

In  Figure  5  the  slit  image  is  shown  sitting  on  top  of  the  rod,  slightly 
overlapping  the  rod  proper.  This  arrangement  implies  the  use  of  the  photo¬ 
graphic  instrument  since  the  spectrum  is  resolved  along  the  slit.  The 
lower  protion  of  the  spectrum  is  then  the  black  body  radiation  from  the 
surface  on  which  is  superimposed  the  local  gaseous  radiation.  The  upper 
part  of  the  spectrum  contains  the  radiation  from  the  boundary  layer  which 
is  in  the  region  of  strong  gradients.  From  the  black  body  radiation  the 
surface  temperature  of  the  rod  may  be  determined  directly  by  noting  the 
wavelength  of  peak  intensity  if  it  is  in  the  sensible  range  of  the  spectro¬ 
graph,  or  by  a  slightly  more  complicated  two  color  analysis  (Appendix  E). 

Since  the  discrete  spectrum  is  from  a  region  of  strong  gradients 
the  radial  variation  of  intensity  must  be  determined  by  a  numerical  solution 
of  Abel’s  integral  equation  (Appendix  F).  Once  the  radial  grandients  have 
been  determined  the  specie  densities  and  temperatures  may  be  calculated 
as  before. 
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Figure  5:  Ablation  Experiment  Setup 


The  actual  flux  rates  of  electrons  from  the  surface  must  be  deter¬ 
mined  from  an  a  priori  estimate  of  local  velocities  since  under  the  con¬ 
ditions  of  the  experiment  the  boundary  layer  velocity  distributions  cannot 
be  measured.  The  thermal  thickness  can  be  measured,  however,  which 
is  of  considerable  help.  The  basic  boundary  layer  equations  to  be  used 
are  given,  in  Appendix  G. 


Since  the  photographic  instrument  does  not  resolve  time,  the  light 
beam  from  the  gas  will  be  split,  half  going  to  the  photographic  instrument 
and  half  to  the  recording  monochromator  which  will  be  set  to  cover  a  line 
of  particular  interest,  e.  g.  ,  a  str|)ng  Cs  or  Cs  line.  In  this  way  the  re¬ 
lative  evolution  rate  may  be  traced  and  the  photographic  data  taken  at  the 
most  opportune  moment.  This  is  important  since  the  samples  have  a 
surface  layer  of  pure  seeding  material  which  will  flash  at  the  start  of  a 
run.  After  the  flash  a  steady  evolution  of  material  occurs  due  to  ablation 
and  it  is  during  this  latter  period  that  the  quantitative  data  must  be  taken. 


The  axial  position  (x)  of  the  measurements  is  yet  to  be  determined 
precisely  and,  since  there  are  fairly  strong  axial  temperature  gradients 
in  the  carbon  rod,  the  measurement  stations  must  be  chosen  with  care.  The 
three  most  likely  places  are  cone  midpoint,  just  aft  of  the  cone  shoulder, 
and  a  far  downstream  point  where  the  rod  is  relatively  cool.  A  standard 
set  of  locations  will  be  selected,  however,  in  order  to  compare  the  data 
empirically. 

The  ablation  samples  to  be  tested  initially  are  1/2”  diameter  carbon 
rods  which  have  been  vacuum  impregnated  with  materials  or  compounds  of 
materials  which  are  known  to  have  a  low  ionization  potential.  These  are 
aluminum,  sodium,  potassium  nitrate  and  cesium  nitrate.  The  reasons 
for  selecting  the  salts  of  cesium  and  potassium  are  to  minimize  the  handling 
hazards  and  to  eliminate  the  formation  of  alkali  carbides  which  have  little 
o  r  no  mechanical  strength. 


Since  the  total  amount  of  impregnated  material  is  known,  a  rough 
estimate  of  the  material  lost  will  be  obtained  by  weighing  the  sample  after 
a  run,  assuming  that  no  carbon  is  lost.  A  quantitative  chemical  analysis 
of  the  ablated  sample  will  give  better  results  and  will  be  performed  if  tech  - 
niques  are  available.  (American  Metal  Products  made  the  original  samples 
and  will  be  sent  the  used  samples  for  chemical  analysis.  ) 
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Z.  3  Kinetics  Experiment 


The  kinetics  experiment  as  described  below  differs  from  the  pre  - 
ceding  experiments  in  that  the  desired  data  is  quantitative  chemical  rate 
data  rather  than  a  demonstration  of  feasibility.  Also,  the  experimental 
setup  only  grossly  simulates  the  augmenter  device.  Indeed,  the  similarity 
is  accidental  rather  than  nec es sary. 

The  basic  problem  of  a  kinetics  experiment  is  the  temporal  re  - 
solution  of  time  -dependent  events.  The  means  by  which  time  resolution 
is  provided  depends  largely  on  the  time  scale  of  the  process  being  investi¬ 
gated.  For  the  case  of  reactions  at  high  temperatures,  su^h  as  those  con¬ 
sidered  here,  the  time  scale  is  of  the  order  of  10"^  to  10  seconds.  Such 
short  times  preclude  the  investigation  of  a  homogeneous  reaction  in  a  closed 
vessel  or  slow  flowing  gases  in  ducts.  Residence  times  of  the  order  of 
10  to  100  microseconds  are  conveniently  obtained  in  DeLaval  nozzles. 
Furthermore,  when  the  sample  gas  is  heated  by  an  electric  arc,  the  initial 
conditions  of  temperature,  enthalpy  and  pressure  are  controllable  at  will 
and  are  not  dependent  upon  the  chemical  heat  release.  It  is  for  these  rea¬ 
sons  that  the  supersonic  flow  in  a  nozzle  has  been  chosen  as  the  time  re¬ 
solution  mechanism  for  these  kinetic  experiments.  The  basic  apparatus 
is  shown  in  a  schematic  form  in  Figure  6.  The  general  method  of  operation 
i  s  as  follows: 

Argon  (or  other  gases  if  desired)  is  injected  into  the  arc  chamber 
at  rates  consistent  with  good  operating  conditions  for  the  arc.  The  arc  is 
struck  and  the  gas  is  heated  to  some  stagnation  temperature  equal  to  or 
greater  than  the  desired  value.  Because  of  the  nature  of  arc  processes 
the  resulting  hot  gas  is  in  a  highly  non -equilibrium  state.  The  gas  then 
flows  to  the  settling  chamber  where  the  long  residence  times  allow  equili¬ 
bration  to  occur.  If  necessary,  additional  cold  gas  may  be  injected  here 
to  reduce  the  stagnation  temperature  to  the  desired  value. 

At  this  time  the  chemical  system  under  consideration  is  also  in¬ 
jected,  in  this  case  CsNO^  and  A1  .  The  quantities  are  so  selected  to 
make  the  total  heat  of  the  reaction  insignificant  compared  to  the  enthalphy 
of  the  diluent.  This  allows  one  to  ignore  the  reaction  thermodynamically 
and  compute  local  temperatures  and  pressures  based  upon  the  diluent 
properties  only.  This  latter  is  not  essential  to  the  experiment  but  the 
simplication  in  the  analysis  is  considerable. 
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Figure  6 

Schematic  of  Flow  System 
For  Kinetics  Experiment 


Because  of  the  comp] ty  and  number  of  possible  reactions  which 
Cjan  occur  with  the  reactants  selected,  it  is  desirable  to  sort  out  empirically 
the  dominant  reactions.  The  procedure  for  doing  this  is  to  vary  the  con¬ 
centrations  of  species  independently  of  what  is  obtained  from  the  initial 
reactants.  All  atomic  species  except  Cs  can  be  varied  at  will  by  mixing 
with  the  prime  gas.  The  various  percentages  which  will  be  added  cannot 
be  pr e-estimated  specifically  but  factors  of  2  or  3  are  centainly  reason¬ 
able.  The  system  of  chemical  and  rate  equations  along  with  present  esti¬ 
mates  of  the  rate  constants  are  given  in  Appendix  I. 

The  most  interesting  ones  are,  of  course,  those  which  involve  electrons. 
The  equations  which  do  not  involve  electrons  are  present  because  they  form 
part  of  the  chain  which  supply  atomic  and  molecular  species  appearing  in 
the  electron  reactions.  The  specie  concentration  measurement  will  be 
made  in  the  same  manner  as  described  for  the  first  experiment.  The 
method  of  analysis  for  simple  reactions  is  described  in  Reference  2  and  re¬ 
peated  in  Appendix  H.  A  similar  method  will  be  used  to  obtain,  a  global 
electron  recombination  rate  for  the  more  complicated  set  of  reactions  in¬ 
volved  in  the  present  process.  Because  of  the  complexity,  the  results  will 
most  probably  show  a  non-integral  order  and  will  be  strictly  applicable 
only  to  the  conditions  of  the  test.  The  test  -conditions  of  Reference  2  per¬ 
mitted  the  use  of  glass -walled  nozzles  and  thus  only  a  single  nozzle  was 
required.  For  the  present  tests,  transparent  walls  do  not  appear  practical. 

To  trace  the  composition  as  a  function  of  position  in  the  no  zzle,  a  series 
of  nozzles  will  be  used,  each  having  a  different  residence  time.  Since  the 
exit  velocity  is  only  a  slowly  varying  function  of  pressure  ratio  for  very 
large  ratios,  it  follows  that  the  residence  time  is  predominantly  a  function 
of  nozzle  length.  For  this  reason,  it  appears  simple  to  fix  pressure  ratio 
and  vary  the  residence  time  by  varying  length  only,  such  that  the  area  is  a 

common  function  of  the  non-dimensional  axial  position  (  '  ~  X/L  )  where  L» 

^  O'  o 

is  the  nozzle  length,  i.  e.  ,  a  nozzle  which  '’stretches'*  in  the  x  direction 
only.  Therefore,  the  time  to  attain  any  pressure  and  temperature  in  the 
nozzle  will  scale  directly  as  the  length.  Also,  a  single  equilibrium  and 
frozen  calculation  will  suffice  for  any  given  set  of  initial  conditions.  The 
method  of  analysis  for  the  more  complicated  case  is  presently  being 
developed. 
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APPENDIX  A 
RESOLUTION  CELL 


It  is  desirable  to  define  the  resolution  cell  of  plasma  referred  to 
the  optical  system  used  in  the  experiment.  The  system  is  arranged  as  in 
Figure  A-1,  where  the  single  lens  shown  represents  the  combined  optical 
system. 

Cylindrical  symmetry  is  assumed  for  the  plasma,  and  the  radiation 
density  is  assumed  to  be  a  function  of  the  radial  distance,  r,  from  the 
axis  of  symmetry.  The  analysis  assumes  a  cylindrical  slab  of  plasma 
of  thickness  h,  where  h  is  a  characteristic  dimension  of  the  spectrograph 
slit.  » 


Assuming  a  radiation  density  variation  with  radial  distance  only, 
the  radiation  emitted  by  the  incremental  volume  is  given  by  ‘ 

dl  =  I  (r)  dV 

=  I  (r)  dr  dA 

=  I  (r)  dr  (hr  dO) 

The  incremental  solid  angle  subtended  by  the  radiating  volxime,  per  imit 
slit  area  §A,  is  given  by 


d(xi 


hr  do 


Before  performing  the  indicated  integration  over  0,  the  expression 
for  that  half  of  the  resolution  cell  lying  between  the  slit  image  and  the 
optics  will  be  shown  to  be  the  same  as  that  already  derived.  This  con¬ 
clusion  is  based  on  the  symmetry  of  the  system  and  follows  from  the  fact 
that  light  entering  the  optical  system  is  treated  as  though  it  had  first 
passed  through  the  slit  .  The  foregoing  development  is  reversed, 

but  the  same  incremental  quantities  are  involved. 


% 
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Now,  performing  the  integration 


and 


so  that 


Si 


dl  =  I  (r)  dr  (hrd) 


&-  =  I  (r)  Sa  dr 


Thus,  it  is  seen  that  each  quasi-annular  increment  of  the  resolution 
cell  is  equivalent  to  an  incremental  volume  of  cross-section  equal  to  the 
slit  area,  having  the  same  radial  variation  in  radiation  density. 
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APPENDIX  B 


CALIBRATION  PROCEDURES 
I.  Calibration  of  Working  Lamps 

Tvingsten  lamps  are  available  which  have  been  calibrated  for 
spectral  radiant  intensity  (watts/ ste radian)  at  specified  currents,  by  the 
National  Bureau  of  Standards.  The  tungsten  lamp,  which  is  very  nearly 
a  blackbody  radiator,  operates  at  a  temperature  which  depends  upon  the 
resistance  of  the  ribbon  filament  and  the  current  passing  through  the 
filament.  Due  to  boiling  away  of  the  tungsten  filament  and  subsequent  in- 
crease  in  filament  resistance  and  coating  of  the  lamp  inner  surface,  the 
spectral  radiant  intensities  of  the  standard  lamps  decrease  as  the  lamps 
are  used.  Unfortunately,  the  high  cost  of  National  Bureau  of  Standards 
calibrated  lamps  makes  their  use  as  "working”  standards  undesirable. 
However,  by  using  the  known  characteristics  of  the  N.  B.  S.  calibrated 
lamps,  the  spectral  radiant  intensity  of  any  tungsten  lamp  may  be  de¬ 
termined.  (It  is  important  that  the  lamps  be  of  identical  composition  in 
order  that  the  grey  body  emissivities  of  the  filaments  and  the  transmis¬ 
sion  of  the  lamp  materials  be  identical. ) 

According  to  Planck* s  law  of  radiation,  the  relative  wavelength 
distribution  of  radiant  flux  depends  only  upon  the  temperature  of  the 
emitter  (see  Figure  B-1).  Also,  Wein's  displacement  law  states  that  the 
wavelength  of  the  maximum  radiance  varies  inversely  with  temperature 
(also  Figure  B-1).  Thus,  for  tungsten  filaments  at  the  same  temperature, 
the  relative  radiant  intensity  curves  (normalized  to  a  given  wavelength) 
will  be  identical.  It  remains  then  to  determine  the  constant  of  propor¬ 
tionality  which  exists  between  the  outputs  of  the  NBS  lamp  and  the  working 
lamp. 

Determination  of  Filament  Temperature 


For  a  given  lamp,  the  filament  temperature  is  a  function  of  the 
current,  the  dependency  being  determined  in  the  laboratory  by  finding  the 
location  of  the  radiant  intensity  maxima  at  various  currents  and  applying 
Wein*s  displacement  law.  That  one  has  determined  an  actual  radiant  in¬ 
tensity  maximum  rather  than  the  maximum  sensitivity  point  of  the 
measuring  device  may  be  checked  by  calibrating  the  measuring  device 
(sensitivity  as  a  function  of  wavelength)  using  the  NBS  lamp  and  the  known 
data  for  that  lamp.  Then,  as  shown  by  the  dashed  lines  in  Figure  B-2, 
it  is  a  simple  matter  to  determine  the  desired  operating  current  for  the 
working  lamp. 
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Figure  B-2;  Filament  Characteristics 


Recording  System  Attenuation  Table 


.761+3  8.690+2  3.513+2  1.774+2  8.976+1  3.653+1  1.863+1  9.446+0  3. 845+0  1. 961+0 


Figure  B-4; 


Ebert  Optics 
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Determination  of  Proportionality  Constant 


The  physical  units  of  the  measurements  are  trivial  at  this  stage. 

It  is  merely  important  that  the  operating  and  measurement  conditions  for 
the  NBS  lamp  and  the  working  lamp  be  identical.  For  equal  temperatures, 
the  lamp  output  is  also  a  function  of  source  area,  as  shown  in  Figure  B~3. 
In  this  way,  the  constant  of  proportionality  is  found  between  the  working 
lamp  and  the  NBS  lamp  such  that  the  spectral  radiant  intensity  of  the 
working  lamp  may  be  determined.  In  general,  this  proportionality  will 
vary  somewhat  with  the  wavelength  due  to  dissimilarities  in  the  lamp  ma¬ 
terials.  It  may  or  may  not  be  important  to  take  this  variation  into  account, 
depending  upon  the  accuracy  of  the  measurements  and  the  extent  of  the 
variation. 

II.  J.  A.  0,  5  m.  Ebert  Scanning  Spectrometer 

The  Ebert  Spectrometer  consists  of  an  entrance  slit,  a  collimating 
concave  mirror,  a  blazed  grating,  and  an  exit  slit  (Figure  B-4)<.  For 
calibration,  the  image  of  the  filament  of  the  working  lamp  is  focused  upon 
the  entrance  slit  and  the  exit  slit  is  focused  upon  a  photo-multiplier  tube 
or  any  other  detecting  device.  The  output  of  the  detector  is  amplified  and 
used  to  drive  a  strip  chart  recorder.  Various  attenuation  factors  are 
built  into  the  electronic  recording  system  of  the  spectrometer  in  order  to 
be  able  to  uniformly  distinguish  spectral  lines  which  are  either  very 
bright  or  very  dim.  These  attenuation  factors  may  be  set  in  by  turning 
a  switch. 

Intensity  Calibration 

In  quantitative  spectrographic  work,  it  is  desirable  to  be  able  to 
determine  the  radiant  intensity  of  a  given  spectral  line.  It  therefore  be¬ 
comes  desirable  to  determine  the  response  of  the  spectrometer  to  all 
wavelengths.  This  will  be  a  function  of  the  spectral  reflectivities  of  the 
optics,  spectral  sensitivity  of  the  detector,  and  other  factors.  The  work  ¬ 
ing  lamp  is  placed  in  the  position  where  the  gas  to  be  analyzed  will  be 
located,  and  the  image  of  the  filament  is  focused  upon  the  entrance  slit. 

It  now  remains  to  scan  through  the  spectrum  of  the  working  lamp  and 
record  the  photo -multiplier  output.  The  procedure  may  require  setting 
in  some  attenuation  with  the  attenuation  switch.  Using  the  following 
table,  the  output  is  divided  by  the  attenuation  fa  ctor  and  then  compared 
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with  the  known  spectral  radiance  of  the  working  lanmp.  Dividing  each  out¬ 
put  signal  by  the  known  radiance  at  the  given  wavelength,  the  chart  de¬ 
flection  per  watt/steradian  is  determined  as  a  function  of  wavelength. 

Notice  that  the  lamp  must  be  run  at  the  current  defined  by  Fig¬ 
ure  B-2,  and  that  a  small  error  in  this  current  will  be  doubled  in  the 
output.  Hence  it  is  necessary  that  the  lamp  current  be  monitored  ac¬ 
curately. 

Wavelength  Calibration 


The  Ebert  is  calibrated  for  wavelength  by  viewing  the  spectrum 
of  a  standard  source,  such  as  an  iron  arc.  The  iron  arc  has  over  300 
lines  in  the  wavelength  region  from  2400  %to  6700  Other  stand¬ 

ards  are  mercury,  neon,  and  cesium  arcs.  The  wavelengths  are  read 
out  directly  by  a  Veeder  root  counter. 

It  may  be  desirable  to  calibrate  for  wavelength  before  calibrating 
for  intensity,  as  there  may  be  effects  (e.  g.  ,  absorption  lines)  which  in¬ 
fluence  intensity  measurements  at  certain  discrete  wavelengths. 

Ill,  J.  A.  1.5  m,  Wadsworth  Spectrograph 

The  Wadsworth  spectrograph  (see  Figure  B-5)  consists  of  an  en¬ 
trance  slit,  a  concave  collimating  mirror,  a  concave  blazed  grating,  and 
a  curved  strip  of  film.  Calibration  of  the  device  is  complicated  by  the 
non-linearity  of  the  wavelength  dispersion  scale  and  the  non-linear  char¬ 
acteristics  of  the  photographic  emulsion. 

Wavelength  Calibration 

Fixed  to  the  film  guide  is  a  wavelength  reference  marker.  This 
shields  part  of  the  film,  leaving  an  unexposed  reference  mark  on  the  de¬ 
veloped  strip  (Figure  B-6).  Wavelength  calibration  is  accomplished  by 
photographing  the  spectrum  of  a  standard  source,  as  with  the  Ebert 
spectrometer.  The  intensity  variation  is  determined  using  a  micro¬ 
densitometer,  and  the  wavelength  calibration  data  is  recorded  as  ^ 
versus  micro-densitometer  travel  from  the  reference  mark. 
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Figure  B-5:  Wadsworth  Optics 


Figure  B-6; 


Wavelength  reference 
Wavelength  Marker 


Film  strip 


Intensity  Calibration 


Before  it  is  possible  to  discuss  the  intensity  calibration  of  the 
instrument,  it  is  necessary  to  mention  some  of  the  characteristics  of 
photographic  emulsions  since  these  are  used  for  recording  the 

spectral  data. 


The  photographic  emulsion  does  not  respond  in  a  linear  fashion. 

If  one  plots  the  density  (logiQlo/^)  versus  the  logjo  of  the  exposure  (in¬ 
tensity  X  time),  a  characteristic  curve  is  obtained  such  as  that  shown  in 
Figure  B-7.  For  the  linear  portion  of  the  curve,  the  density  may  be 
given  as  a  function  of  the  exposure  raised  to  the  /  power.  All  photo¬ 
graphic  emulsions  exhibit  dependence  upon  the  wavelengths  with  which 
they  are  illuminated.  Most  have  a  variation  in  with  \  for  a  given 
photographic  process  (Figure  B-8a).  All  of  them  have  a  strong  dependency 
with  regard  to  sensitivity  as  a  function  of  X  (Figure  B-8b).  In  any  case, 
the  selection  of  an  emulsion  depends  upon  the  wavelength  region  to  be  in¬ 
vestigated,  and  it  may  be  necessary  to  use  several  emulsions  to  cover  the 
desired  wavelength  range.  Thus,  the  selection  of  suitable  photographic 
emulsion  is  the  first  step  in  the  intensity  calibration  procedure.  By 
careful  selection  of  a  development  procedure,  it  is  possible  to  make  'y 
equal  to  one,  in  which  case  the  density  is  directly  proportional  to  the 
exposure.  Once  the  development  process  has  been  fixed,  one  finds  it 
desirable  to  operate  in  the  range  of  exposures  in  which  the  characteristic 
curve  is  linear.  Thus,  the  determination  of  a  standard  development 
technique  for  the  chosen  emulsion  is  the  second  step  in  the  intensity  cali¬ 
bration  procedure. 


A  final  consideration  in  the  discussion  of  emulsions  is  the  recip¬ 
rocity  law  of  Bunsen  and  Roscoe,  which  states  that  the  exposure  is  equal 
to  the  product  of  the  intensity  and  the  exposure  time.  It  has  been  shown 
that  equal  exposures  do  not  produce  equal  effects  on  the  film  (reciprocity- 
law  failure).  Furthermore,  when  (as  with  the  Wadsworth)  the  exposure 
along  the  slit  is  varied  by  spinning  a  stepped  sector  over  the  entrance  slit 
(Figure  B-9),  the  intermittent  exposures  do  not  produce  the  same  effect 
as  equivalent  continuous  exposures  except  above  a  critical  frequency 
(intermittency  effect).  For  most  films  a  safe  operating  frequency  is 
about  1600  rpm,  although  with  some  emulsions  (and  this  is  important) 
there  is  always  an  intermittency  effect  with  the  sector  arrangement 
It  is  important  that  the  illumination  is  constant  over  the  slit. 
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Hence,  the  ribbon  filament  is  especially  desirable  in  calibration  work 
since  it  gives  relatively  viniform  illumination  over  the  entire  entrance 
slit. 


Once  an  emulsion  is  chosen,  it  is  placed  in  the  spectrograph 
and  the  working  lamp  spectrum  is  photographed.  Micro-densitometer 
traces  are  made  of  the  resulting  negatives,  and  a  similar  procedure  is 
followed  to  determine  the  wavelength  response  as  for  the  Ebert  spectrom¬ 
eter.  It  is  important  here  to  try  several  exposure  conditions  to  determine 
if  one  is  operating  in  the  linear  region  of  the  film. 
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APPENDIX  C 


MEASUREMENT  OF  ION  DENSITIES  USING 
THE  HOLTZMARK  THEORY  OF  LINE  BROADENING 


Theory 

It  is  shown  by  Holtzmark  )  that  the  average  electric  field  at 
a  radiating  atom  in  a  plasma  due  to  the  ions  alone  is  (for  singly  charged 
ions) : 

2/3 

Eq  =  2.6len^  (C-1) 


where  e  =  electronic  charge  and  n^  =  ion  density.  The  statistical  dis¬ 
tribution  of  the  electric  field  (isotropic)  is  also  shown  to  be 


W 


(^)  d^  =  1.496^ 


-5/2 


-3/2 


(1  +  5.  106/?  + 
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43,;?"^+. 


.  )  d/-?  (C-2) 


dy5  1^^  large  j 

E  r  ' 

where  ^  =  "g —  =  (  — ) 


(C-2a) 


(C-3) 


and  r  is  the  distance  between  the  ion  and  radiator. 

Note  that  large  imples  snaall  >  (i.  e.  ,  close  encounters)  and  the  latter 
form  (C-2a)  applies  for  the  large  E  (i.  e.  ,  in  the  wings  of  the  spectral 
line) . 


The  first  order  Stark  effect  gives  the  displacement  of  the  per¬ 
turbed  line  from  the  original  undisturbed  line  center  as  being  proportional 
to  the  electric  field  strength: 

A  ^ 
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where  A  X  the  line  shift  (in  angstroms)  and  is  the  proportionality 
constant  for  the  particular  component.  Thus 


(C-5) 


If  the  unperturbed  spectral  line  consists  of  several  components  whose 
intensity  and  Stark  coefficients  differ,  then  the  observed  line  will  consist 
of  the  sum  of  terms 


where  J  — J  d  ^  represents  the  line  intensity  in  the  interval 

d  (A  A)  about  A  A  and  is  the  transition  probability  of  the  component 
of  the  line  in  question. 

.  aa 

Equation  (C-6)  may  be  written  in  terms  of  the  parameter 


V  CL 

J  (a)  da  =  /  w  (  -^)  da 


(C-7) 


k 


J  (CL)  is  readily  computed  using  the  known  values  of  Aj^  ,  ,  and 

Equation  (C-2). 


The  parameter  Cl  has  been  defined  as: 


AX.  °  , 

-r=. —  in  A  kilovolt 


(C-8) 


where  A/'  is  the  displacement  from  the  line  center  and  Eq  is  the  normalized 
electric  field  strength  which  is  related  to  the  electron  density  by: 

3/2 

(C-9) 


n;  - 


2.  61 


e  -I 
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where  e  =  electronic  charge  =  4.8025  x  10 


esu 


-10 


Eq  =  normalized  electric  field 
=  esu  volts  =  300  volts 


=  0.3  kilovolts 


from  the  definition  of  U.  it  is  seen  that* 


Qj  -0^2 

~a~i 


=  C  (C-IO) 


The  theoretical  profiles  have  been  replotted  in  Figures  C-1,  C-2,  and 
C-3  to  give 


If  one  plots 


a 


1 


theo 


vs.  ClLj  for  C  =•  0.  5 


a. 


R 


AA, 


R 
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exp 


R 
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^R 


a2 


(C-ll) 


VS.  A  X  2  for  - r -  =  C  =  0.  5,  as  in  Figures  C-4  and  C-5,  the  proper 

1 

value  of  and  is  obtained  immediately  from  the  second  equality 

in  (C-ll). 


Application  of  (C-1)  and  (C-2)  yields  the  required  electron  concentration. 

Because  of  the  limitations  of  the  theory,  the  method  is  not  applicable 
near  the  center  of  the  line  and  no  attempt  should  be  made  to  fit  the  center. 

A  little  experience  will  give  the  necessary  feel  for  the  best  fit. 

Because  of  the  difficulty  of  reading  the  experimental  data,  Figure 
C-6,  at  the  wings  of  the  line,  the  profile  may  be  plotted  on  log-log  paper 
and  extrapolated  with  a  straight  line.  The  slope  will  usually  be  found  to 
be  -5/2  as  in  Figure  C-7  and  the  required  data  read  off. 
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TO  X  to  TO  THE  CM.  309-14 
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For  this  analysis  the  experimental  data  need  not  be  corrected  for 
phototube  response  since  the  response  will  not  vary  appreciably  within  the 
wavelength  span  of  a  single  line  and  only  relative  values  of  amplitudes  are 
required.  In  practice,  the  result  from  a  single  line  should  not  be  con¬ 
sidered  conclusive.  Comparison  with  the  results  of  other  lines  should  be 
made.  In  general  the  results  may  differ  by  a  factor  of  2  or  so.  The  final 
result  depends  upon  the  quality  of  the  data  and  the  theoretical  fit. 

Note  that  while  the  theory  apparently  indicates  a  shift  in  one  di¬ 
rection  only,  the  quantum  mechanics  of  the  Stark  effect  shows  a  symmetrical 
splitting,  and  for  linear  Stark  effect  the  center  of  the  line  is  not  displaced. 

If  desired,  both  sides  of  the  measured  line  may  be  used  for  analysis  and 
the  results  averaged  since  it  is  sometimes  difficult  to  read  the  precise 
center  from  the  data.  Occasionally,  for  broad  lines,  the  background  or 
other  close  lines  will  dictate  the  most  convenient  half  of  the  line  to  be  used. 
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APPENDIX  D 


ALTERNATE  METHOD  OF  ION  DENSITY  MEASUREMENT 


The  results  of  a  more  rigorous  derivation  of  the  broadening  effect 
of  ions  and  electrons  is  given  by  Griem^^"^^  in  which  the  half  width  at  half 
height  of  several  lines  in  the  visible  cesium  spectrum  are  shown  as  a  linear 
function  of  electron  density  with  a  slight  dependence  on  temperature.  The 
theoretical  details  for  Griem's  calculations  are  not  available  and  hence  the 
form  of  the  terrperature  dependence  is  not  known.  In  view  of  the  small 
variations  and  other  uncertainties  it  will  be  assumed  that  linear  interpolation 
is  sufficient.  This  should  certainly  lead  to  errors  no  larger  than  the  rough 
estimate  of  10%  to  Z0%  of  the  ion  broadening  effect. 

The  numerical  data  is  given  in  Appendix  J. 

The  data  reduction  procedure  is  extremely  simple: 


1.  Determine  the  temperature  by  the  methods  of  Appendix  E. 

Z.  Measure  the  half  width  at  half  height  for  the  prominent  Cs  lines. 
(See  Figure  D-1) 


3.  Interpolate  between  tabxilated  values  to  obtain  theoretical 
broadening  at  the  measured  temperature. 


4.  Scale  the  electron  density  directly  with  line  width  with  a  0.8  to 
0.  9  reduction  due  to  the  ion  broadening  effect. 

to  o 


CO  o  _ 

2500 


5000  theoretical  broadening 

values  at  the  reference  electron  density  (10^^/cc)  and  temperature,  then 


(f-X 


2500 


C-)  -(-) 

^^“^  5000  ^2/ 


2500 


/T-2500\ 
\  2500  y 


is  the  reference  broadening  at  the  measured  temperature  T  (°K). 


The  measured  electron  density  is  then: 
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I 

I 

1 


(f), 

where  ^  is  a  parameter  accounting  for  the  ion  broadening  and  is  typically 
0.  8  to  0.  9. 
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APPENDIX  E 


OPTICAL  METHODS  OF  TEMPERATURE  DETERMINATION 

Temperature  profiles  must  be  determined  in  two  specific  regions: 
in  all  phases  of  the  program,  gas  temperature  levels  are  required;  in  the 
ablation  experiment,  the  surface  temperature  of  the  test  model  must  also 
be  obtained. 


Black  Body  Radiation 


The  radiation  from  the  surface  of  the  ablation  test  model  is  con¬ 
sidered  to  be  characteristic  of  a  black  body  at  the  surface  temperature  of 
the  model.  Thus,  the  intensity  is  given  by: 


(E-l) 


where  I  is  the  radiant  intensity  in  units  of  power  per  unit  area, 


K  is  constant  depending  on  the  optics  of  the  system  and  the  emissivity  of 
the  body  under  analysis,  cj^,  C2  Sire  the  first  and  second  radiation  constants, 
respectively,  X  Is  the  specific  wavelength  at  which  the  intensity  measure¬ 
ment  is  made,  and  T  is  the  temperature  of  the  radiating  surface. 


Consideration  of  the  a  priori  limits  imposed  on  the  experiment 
allows  a  further  simplification  of  Equation  (E-l).  Thus,  the  surface 
temperature  of  the  test  model  can  be  expected  not  to  exceed  4000°K 
(melting  point  of  graphite  3900^K),  and  the  response  of  the  phototube  in 
the  optical  system  limits  the  useable  wavelength  region  to  4000  K  to 
9000  nominally.  Substitution  of  these  maximum  values  (i.e.  X  ^ 

36  cm°K)  into  Equation  (E-l)  shows  that,  to  within  2%. 

Ix  =  K  ^  exp  (-  ^  )  (E-la) 

^  x^ 

The  evaluation  of  local  surface  temperature  proceeds  from  the  relative 
intensities,  measured  at  two  different  wavelengths,  radiating  from  the 
same  element  of  body  surface  area.  The  ratio  of  two  intensity  measure¬ 
ments  is  given  by: 
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5 


(E^2) 


XI 


^  2 


The  temperature  is  then  given  by: 


T 


^2 


\  -  ^^2 

1 

j^l  >  2 

J 

+  5  In 


-1 


It  is  seen  from  Equation  (E-3)  that,  as  approaches  Xp  T  becomes 
the  ratio  of  two  numbers  close  to  zero.  Thus,  to  ensure  better  accuracy, 
the  comparison  of  measured  intensities  should  be  carried  out  for  wave¬ 
lengths  as  far  apart  as  possible.  Also,  the  ratios  of  several  measure¬ 
ments  should  be  used,  and  a  mean  value  of  the  computed  temperature 
determined. 


The  value  of  cz  required  is  dependent  on  the  units  used  for  X  •  For  \  in 
Angstrom  units, 


cy  =  ^  =  1.4387  x  10®  A  -  °K 
2  k 


The  expected  error  in  a  single  temperature  value  cannot  be  estimated  a 
priori,  due  to  the  complicated  dependence  of  T  on  both  the  measured  in¬ 
tensity  and  the  wavelength  at  which  the  intensity  measurement  is  made. 
Differentiating  Equation  (E-3),  a  formula  for  the  relative  error  in  T,  as 
a  function  of  the  uncertainty  in  the  measured  intensity  is : 


At 

T 


-1 

Ai  1  ^ 

^1  h 

L  J 

(E-4) 


Thus,  the  experimental  data  will  allow  an  estimation  of  the  theoretical 
error. 


Gas  Phase  -  Line  Spectra 

The  spectral  radiation  from  several  gas  constituents  will  be  used  to 
determine  the  temperature  of  the  flow  both  around  the  ablation  test  model 
and  from  CSNO3:  A1  rocket  motor.  The  emitting  species  in  the  rocket 
motor  exhaust  will  be  cesium  and  aluminum.  In  the  ablation  tests,  sodium 
and  potassium  will  be  used  as  well.  In  addition,  it  may  be  necessary  to  seed 
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the  argon  of  the  plasma  jet  with  hydrogen  in  order  to  obtain  emission  in 
the  proper  spectral  region.  The  theoretical  basis  for  data  reduction,  as 
well  as  the  requisite  data  for  the  radiating  species, are  given  below. 


I 

I 

1 


The  measured  intensity  of  a  spectral  line  is  given  by; 


li  =  KAi  ^  Ni 


(E-5) 


where  K  is  a  constant  of  the  optical  system, 

is  the  transition  probability  of  the  transition  under  consideration, 
h  is  Planck's  constant, 


c  is  the  speed  of  light  in  vacuo, 

is  the  wavelength  of  emitted  radiation,  and 

is  the  number  density  of  particles  of  the  radiating  specie  in 
the  upper  state  of  the  transition  in  question. 

is  assumed  to  obey  the  Maxwell- Boltzmann  distribution  law,  i.e. 


Ni  =  N 


J  J 


(E-6) 


Ei/kT 


where  N  is  the  total  number  density  of  the  radiating  specie, 

gi  is  the  degeneracy  of  the  upper  state  of  the  transition, 
Ej^  is  the  energy  of  this  upper  state, 
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k  is  the  Boltzmann  constant, 

T  is  the  local  temperature,  and 

Q,  as  defined  in  Equation  6,  is  the  partition  function  of  the  radiating 
specie. 

By  direct  substitution  of  Equation  (E-6)  into  Equation  (E-5)  ,  we  have 

he  N  -E/kT  (E-7) 

X  'q 

or 

=  ln(Khc  -^1-  f-  i  (E-7a) 

IX 

Thus,  it  is  seen  that  if  we  plot  1  n  ("^g) -versus  (E/k)  (i.e.  ,  the  energy 
level  expressed  as  an  equivalent  temperature)  the  magnitude  of  the  slope 
of  the  resulting  straight  line  plot  will  be  the  reciprocal  of  the  desired 
temperature  value,  algebraically: 


T 


d 

d(E/k) 


1 


(In 


(E-8) 


The  slope  of  the  line  may  be  approximately  determined  by  an  empirical 
'‘best  fit"  or  by  a  least  squares  fit  on  a  computer.  A  check  on  the  results 
may  be  obtained  by  using  a  relative  intensity  method  similar  to  that  described 
for  the  black-body  data  reduction. 

The  relative  error  in  the  temperature  due  to  experimental  uncertainty 
in  the  intensity  values  may  be  obtained  from  Equation  (E-7)  by  differentation. 
Thus 


I  = 

Ai 


he  N 
T  "Q 


he  N 

=  X  Q 


■E/kT 


E  At 


-E/kT 

e 


(E-7) 

(E-7b) 
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and 


Ai  ^  At 

I  kT  T 

or 

At  ^  ^ 

T  El. 


(E-9) 


(E-9a) 


Therefore,  the  variation  in  T  may  be  determined  from  the  experimental 
data  for  each  specific  line. 

The  various  constants  required  in  the  above  equations  are  presented 
in  Appendix  J. 
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APPENDIX  F 


DETERMINATION  OF  RADIATING  SPECIE  CONCENTRATIONS  BY  A 
NUMERICAL  SOLUTION  OF  ABEL'S  INTEGRAL  EQUATION 

The  method  outlined  below  was  presented  by  Pearce  '  ~  ^ 

The  data  to  be  reduced  will  be  in  the  form  of  photographic  plates  from  which 
relative  radiant  intensity  may  be  determined. 


It  is  assumed  that  the  plasma  is  cylindrically  symmetric, 
and  that  the  concentration  of  radiating  species  varies  as  a  radial  step- 
function.  (See  Figure  F-1).  The  photographic  raw  data  is  then  divided 
into  an  arbitrary  number  of  equally  spaced  segments  on  this  basis. 


The  analysis  begins  with  the  last  segment  of  the  photographic 
data.  The  intensity,  is  known  and  we  wish  to  find  the  radiating  specie 
concentration.  Thus,  from  the  figure, 


I 

n 


A  A_ 


n,  n 


N 

n 


Aa  is  found  in  the  referenced  work  in  tabular  form  for  n 
less  than  or  equaf  to  25. 


Proceeding, 


n-1  -'n 

where  the  unknown  is  now  N 


I  1  =  N  (Aa  1  )  +  N 

"  n-i.n'  ^ 


n-1  ^  ^  n-1^ 


n-r 


In  general,  is  found  from 
n 


L  = 
J 


(j  -  n,  n-  1,  .  .  .  2  ,  1) 


The  analysis  thus  produces  one  annular  concentration  value  at  a  time, 
beginning  at  the  outside  of  the  plasma,  and  working  into  the  central  core. 
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APPENDIX  G 


BOUNDARY  LAYER  VELOCITY  PROFILE 


It  is  the  purpose  of  this  appendix  to  describe  an  approximate  method 
for  deducing  the  boundary  layer  velocity  profile  from  corresponding  thermal 
boundary  layer  data.  The  latter  is  to  be  obtained  from  plasma-jet  tests  of 
either  blunted  or  conically  tipped  cylinders  using  argon  gas  as  the  basic 
working  medium. 

It  is  anticipated  that  the  plasma-jet  test  facility  will  be  such  that  the 
argon  gas  will  be  chemically  frozen  after  expansion  from  its  initial  reservoir 
(stagnation)  conditions  prior  to  flow  about  the  test  body.  As  a  result,  the 
associated  aerodynamic  flow  processes  are  assumed  to  be  calorically  per¬ 
fect. 


The  final  choice  of  the  body  shape  will  most  likely  be  dictated  by 
the  existing  condition  of  the  test  model  material,  namely  right  cylindrical 
shapes  of  high  slenderness  ratios.  Although  such  bodies  are  difficult  to 
analyze,  little  would  be  gained  from  the  use  of  other  shapes  due  to  lack  of 
simple  analytical  solutions  (or  experimental  data  corresponding  to  antici¬ 
pated  conditions)  to  their  associated  flow  fields.  Further,  if  the  region 
under  examination  lies  several  body  diameters  downstream  of  the  frontal 
nose- cylinder  junction  then  the  inviscid  flow  will  approximate  free  stream 
conditions  regardless  of  the  nose  shape.  Because  of  its  inherent  simplicity, 
such  an  arrangement  will  be  assumed. 

Normally,  simple  solutions  to  the  viscous  boundary  layer  flow  do 
not  exist.  However,  under  certain  restrictive  assumptions,  reasonable 
relationships  may  be  obtained  which  are  directly  applicable  to  the  expressed 
problem.  The  Crocco  analysis  for  compressible  viscous  flow  with  a  Prandtl 
number  of  unity  is  such  a  case  Without  going  into  the  details,  the 

basic  assumptions  are  tabulated  below: 

(a)  No  pressure  gradient  along  the  body  surface 

(b)  No  temperature  gradient  along  the  body  surface 

(c)  Steady-  state  conditions 

(d)  Constant  Prandtl  number. 
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If  the  test  point  is  located  far  downstream  of  the  nose  the  first  two  assump¬ 
tions  will  be  satisfied.  The  validity  of  the  third  assumption  depends  on 
the  length  of  the  test  and  model  size.  The  final  assumption  is  rigorously 
invalid  but  appears  to  approximate  actual  conditions  as  the  Prandtl  number 
ranges  from  roughly  0.  6  at  0°C  to  1.  1  at  1200°C.  The  analysis  yields 
the  following  simple  relationship  between  velocity  and  enthalpy: 

u  ^  u 

hB  =  hw  -  (hw  -  hA)  - "  ^e)  (tJ^) 

For  constant  specific  heats  this  reduces  to: 

/  u  \  ^ 

(^A  "  '^W)  E  )  "  Cue 


u 


Solving  for 


E 


u 


Ui 


-  T-yy)  +  VcTa  -  -4  (Ta  -  ~  '^w) 

2(Ta-  Te) 


where 


Ue  = 

U  = 

Te  = 
Tw  = 
Tb  = 


boundary  layer  edge  velocity  (in  this  case  Ug  is  assumed 
equal  to  the  free  stream  velocity) 

local  boundary  layer  velocity 

boundary  layer  edge  temperature 

wall  temperature 

local  boundary  layer  temperature 

adiabatic  temperature  (assumed  equal  to  reservoir 
temperature) 
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Thus  from  a  knowledge  of  the  test  model  wall  temperature,  boundary 
layer  edge  temperature  (approximated  by  free  stream  value),  the 
reservoir  temperature  and  local  temperature,  an  estimate  of  the  local 
velocity  (normalized  to  the  edge  value  which  is  taken  as  the  free  stream 
value)  can  be  made. 
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APPENDIX  H 


CHEMICAL  REACTIONS,  RATE  EQUATIONS,  AND  REACTION  RATE 
CONSTANTS  FOR  KINETICS  STUDY 

The  chemical  reactions  given  below  are  considered  repre¬ 
sentative  of  the  reacting  system  to  be  considered  in  the  kinetics  experi¬ 
ment.  The  data  have  been  previously  published  H-2)^ 

presented  here  for  convenience.  The  first  set  of  equations  defines  the 
reactions  of  the  system.  The  second  set  expresses  the  rate  of  formation 
of  the  various  species.  The  specific  reaction  rate  constants,  gathered 
from  various  sources  '  ^  presented  in  tabular  form. 


A)  Reactions 

+ 

i.  C  s  +  e  +  M 

k2 


Cs  +  M 


2 .  C  "i"  6 


+  3 

®  ( — 


2Cs 


3.  +  e  +  M 


k6 


O^'  +  M 


+  -4 

4.  +  Cs  +  e  ^  2N  +  Cs 

^8 


5.  2N  +  NO 


29 

klO 


+  NO  +  e 


6.  20  +  M 


7.  2N  +  M 


11 

^12 

^13 

kl4 


O  +  M 
2 


+  M 
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8.  N  +  O  +  M 


NO  +  M 


9.  N  +  NO 


10.  N  +  O^ 


N2  +  O 


NO  +  O 


B)  Rate  Equations 

1-  =  -ki  (Cs^)(e)(M)  +k2(Cs)(M) 

-k3  (Cs2)(e)  +  k4(Cs)^ 

-k5  (02)(e)(M)  +k^(02’)(M) 

-ky  (N2){e)(Cs)  +  k8(N)^  (Cs) 
-1^10  (N2)(e){NO^)  +  k9  (N)^(NO) 

2.  =  -k,  (Cs)(e)(M)  +  k2  (Cs)(M) 

dt 

-ky  (N2)(Cs)(e)  +  kg  (N)^  (Cs) 

3.  I  +k  (Cs)(e){M)  -k2  (Cs){M) 

+  2k3  (Cs2)(e)  -  2k4  {Cs)^ 

+  k7  (N2)(Cs)(e)  -  kg  (N)2(Cs) 

+  + 

4.  ^(^^2)  =  .k3  {Cs2){e)  +  k4{Cs)^ 

dt 

.  d(N2)  -  +  7 

5;  -k^  (N2)(Cs)(e)  +kg  (Nr  (Cs) 

+  kg  (N)^  (NO)  -  k^Q  (N^)  {NO)(e) 

+  ki3  (N)^{M)  -  k^4  (N2)(M) 


k^^  (N)(NO)  -  ki8  (N2)(0) 
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6. 


d(N) 

dt 


7. 


8. 


9. 


10. 


:  2k.^(N2)(ct)(e)  -  2kg  (N)^  (Cs) 

+ 

-2kg(N)2  (NO)  +  2k^^(N2)(NO)(e) 
-2k^^  (N)^  (M)  +  2k^^  (N2)(M) 
-k^^(N)(0)(M)  +  k^^  (NO)(M) 
-k^^(N)(NO)  +  k^g  (N^)(0) 

-k^^  (NKO^)  +  k^o  (N0){0) 

^9-1  .-kQ(N)^(NO)  +  kio(N2)(NO)(e) 

+  kj^5  (N)(0)(M)  -k^^  (NO)(M) 
-k^.^  (N)(NO)  +  k^g  (N^)(0) 

+  k^g(N){02)  -  k^Q  (N0)(0) 

+ 

d(NO)  .  (N)2(N0)  -k  (N  )(NO)(e) 

dt  "  9  10  2 

^^^2^  =  -k^  (0^)(e)(M)  +  kg(02')(M) 

dt 

+  k^^  (O)^(M)  -  k^2 
-k  (O  )(N)  +  k2o(NO){0) 

i  7  w 

*^^^2  ^  -  k.  (O  ){e){M)  -  k^(0  ■)(M) 

2 
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-kj5  (N)(0){M)  +  k^^(NO)(M) 


+k  (N)(NO)  -k  (N  )(0) 
18  2 

+  (N)(0^)  -  k^^  (N0)(0) 


i 

12.  - 
dt 

The  reaction 

ki 

0 

rate  constants  are  of  the  form 

-bi  -Ci/T 

r  T  e  ,  where 

a.  bi 

1 

5. 

X 

10-26 

1.  0 

0 

2 

1.  25 

X 

10'^° 

0 

45200 

3 

1.  6 

X 

10-2 

1.  5 

0 

4 

1.  18 

X 

10-10 

1.  5 

33000 

5 

1.  2 

X 

10-28 

0.  5 

0 

6 

1.  58 

X 

10"^^ 

-1.0 

5340 

7 

4.  58 

X 

j^O-22 

1.  5 

68070 

8 

1. 

X 

10'^^ 

0 

0 

9 

1. 

X 

10-29 

0 

0 

10 

1.  08 

X 

10"^9 

1.  5 

5880 

1 1 

1.  30 

X 

10-^9 

1.  0 

0 

12 

9.  07 

X 

lo-'^ 

1.  5 

59480 

13 

9.  22 

X 

lO'^”^ 

1.  5 

0 

14 

1.  02 

X 

10-i 

1.  5 

113250 

15 

3.  97 

X 

10"2® 

].  5 

0 

16 

9.  25 

X 

10"^ 

1.  5 

75500 

17 

1.  72 

X 

10"^^ 

0 

0 

18 

8.  3 

X 

10"^^ 

0 

37900 

19 

1.  68 

X 

10"^2 

-0.  5 

3120 

20 

5.  07 

X 

10-^9 

0 
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APPENDIX  I 


ANALYSIS  OF  RADIATION  FOR  CHEMICAL  RATE  DATA 


Since  the  major  object  of  the  planned  experiments  is  the  evaluation 
of  the  magnitude  of  the  chemical  kinetics  it  is  necessary  that  the  method 
of  data  analysis  for  chemical  rates  be  described  in  some  detail. 

In  what  follows  it  is  assumed  that  the  concentrations  of  the  important 
species  can  be  measured  to  the  desired  accuracy.  The  accuracy  actually 
obtainable  is  discussed  elsewhere. 

The  major  assumptions  are: 

1.  The  thermodynamic  state-of-the«-nozzle  flow  is  known  a  priori 
implying : 

a.  The  flow  is  near  frozen  or  near  equilibrium  and  the  deviations 
from  either  of  these  conditions  do  not  significantly  alter  the 
thermodynamic  state,  or, 

b.  The  reactants  are  present  in  sufficient  dilution  such  that 
they  do  not  effect  the  thermodynamic  parameters  regard¬ 
less  of  the  extent  of  the  reaction  and  an  ideal  gas  law  may 
be  used. 

Z.  The  reactions  given  in  Appendix  H  adequately  describe  the 
system. 

3.  The  isothermal  rate  expressions  are  sufficiently  accurate  for 
the  generation  terms. 

The  basic  equation  is  the  specie  continuity  equation,  which  can  be 
written  in  the  absence  of  diffusion  as: 

A.  (  in  T)  =  S£L 

dt  gt  (I-l) 


dCi  _  Cj 

dt  y-1 


6  3 


where 


number 


Ci 

T 

t 


=  concentration  of  the  specie, 


cc 


7 

Sci 


local  temperature,  K 
time,  seconds 
isentropic  exponent 

local  rate  of  production  of  i^-*^  specie  due  to  reaction 


=  0  for  frozen  flow 


For  equilibrium  flow 


Sci 


/  0  (as  is  the  case  for  isothermal 


Tt  ■  ■ 

equilibrium).  The  production  of  species  equilibrium  nozzle  flow  is  due 
entirely  to  the  temperature  dependence  of  the  equilibrium  constant,  i.e.  , 
for  the  reaction. 


+  A' 


the  rate  of  depletion  of  electrons  due  to  equilibrium  expansion  is 
SlnCe"  __  1  d  (InKe)  dT 


dt 


dt 


(1-2) 


03) 


It  should  be  noted  that  similar  terms  should  be  included  in  the  near  frozen 
calculations.  However,  under  the  near  frozen  assumption,  it  is  implied 
that  the  rate  of  change  of  concentration  due  to  reaction  is  small  compared 
to  the  change  due  to  expansions.  This  allows  the  use  of  isothermal  form 
of  the  rate  equations.  In  any  case,  the  accuracy  of  any  kinetics  experiment 
does  not  warrant  consideration  of  what  are  essentially  second  order  terms 
in  the  Taylor  expansion  of  the  rate  of  production  except  near  equilibrium 
where  the  first  order  terms  are  zero. 

Returning  to  Equation  (I-l)  and  assuming  that  the  right  hand  side  is 
adequately  represented  by  the  isothermal  rate  of  production  as  given  in 
Appendix  H  we  may  write  for  the  general  case  a  set  of  1 1  equations  of  the 
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form  given  in  (I- 1)  corresponding  to  the  11  reacting  species.  The  unknowns 
are  the  10  rate  constants  of  the  ten  stoichiometric  chemical  equations.  This 
implies  that  the  eleven  equations  are  not  independent  and  in  fact  they  are  not 
since  by  charge  conservation  we  have 

S(e)  0(Cg'^)  S(Cs2^)  S(NO+)  S(02’)  _  „ 

Thus  any  one  of  the  conservation  equations  corresponding  to  the  terms  in 
(1-4)  may  be  eliminated  from  the  set  of  simultaneous  equations  and  used 
later  as  a  consistency  check. 


The  measured  quantities  in  the  experimental  set  up  are  the  specie 
concentrations  at  several  positions  in  the  nozzle.  In  practice,  this  is  done 
by  making  the  measurements  at  the  exits  of  a  family  of  nozzles  related  as 
shown  in  the  sketch. 


I 


I 


4 


Equation  (I-l)  may  be  transformed  tospace  coordinates  by 


dx 

V 


dt 


dCi  Ci  V  d  (  1  nT)_  8ci 

dx  "  y-  1  dT  “  gt 


(1-5) 
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Assuming  that  v  and  T  are  predetermined  functions  of  x  and  all 
the  Ci  are  determined  experimentally  as  a  function  of  x,  the  left  hand 
side  of  (1-5)  is  determined  for  every  point  of  the  nozzle.  Thus  the  right 
hand  side  of  (1-5)  forms  a  set  of  linear  simultaneous  algebraic  equations 
in  the  unknown  rate  constants  for  every  point  along  the  axis  of  the  nozzle 
and  hence  for  the  range  of  temperatures  encountered  in  the  nozzle.  The 
specific  rate  constants  may  be  solved  for,  plotted  versus  temperature  in 
the  usual  way,  and  analyzed  for  the  exponential  and  pre-exponential  factors 
by  well  established  methods. 

Intrinsically  the  method  of  analysis  is  simple  and  correct.  However, 
there  are  two  difficulties  which  are  likely  to  be  encountered;  one  numerical 
and  one  physical. 

The  numerical  source  of  error  is  in  the  left  hand  side  of  (1-5)  which 
is  essentially  the  difference  of  two  slopes;  i.e.  ,  the  total  rate  of  change  of 
concentration  and  the  rate  of  change  due  to  the  expansion  process.  Near 
frozen  flow  these  slopes  are  nearly  equal,  leading  to  the  small  difference 
of  large  numbers. 

Related  to  this,  but  inherent  in  the  diagnostic  instrumentation,  is 
the  accuracy  with  which  the  concentrations  can  be  determined. 

Depending  upon  the  specific  results  of  the  experiment  it  may  be 
necessary  to  use  a  finite  difference  form  of  (I-?'  over  a  substantial  nozzle 
segment,  in  which  case  a  detailed  correlation  with  temperature  is  no  longer 
possible.  For  recombination  reactions  this  defect  is  slight  since  the  tem¬ 
perature  dependence  is  usually  a  small  power  of  T  rather  than  exponential. 
It  is,  however,  just  the  recombination  reactions  which  are  most  interesting 
for  the  present  problem. 


66 


APPENDIX  J 


PHYSICAL  DATA 

It  is  the  purpose  of  this  appendix  to  set  down  in  one  place  as  much 
of  the  required  physical  data  as  possible,  in  order  to  facilitate  the  reduc¬ 
tion  of  experimental  data.  The  following  material  has  been  obtained  from 
many  sources  J-2,  J-3,  and  includes  constants,  conversion 

factors,  and  spectroscopic  data  for  the  species  expected  to  be  present  in 
the  gases  analyzed. 

1 )  Physical  Constants 

-27 

h  =  Planck*s  constant  =  6.6251  x  10"  erg-sec. 

c  =  Light  velocity  =  2.9978  x  10^^  cm/ sec 

k  =  Boltzmann* s  constant  =  1.  3805  x  10“^^  erg/°K 

2 

=  First  radiation  constant  =  Z'jfhc 
=  3.  7402  X  10"^  erg-cm^/sec 

c^  =  Second  radiation  constant  =  hc/k 
=  1.4388  cm  - 

e  =  Electronic  charge  =  4.  8025  x  10"^®  esu 

23 

Aq  =  Avogadro's  number  =  6.0228  x  10  particles/mole 

2)  Energy  Conversion  Factors 


erg 

cm"  ^ 

electron  volt 

erg 

1 

5.  0358x  10^^ 

6.  2421 X 10^^ 

cm"  ^ 

1. 98258  X  10^^ 

1 

1. 2395x10“^ 

electron  volt 

1. 60203  X  10"^^ 

8067. 5 

1 

67 


3)  Spectroscopic  Data 
A)  Atomic  Specie 

i)  Nitrogen  and  Oxygen 

The  following  spectral  lines  due  to  atomic  oxygen  and  nitrogen 
may  be  observed: 


Oxygen 

Nitrogen 

o 

o 

4370  A 

4110  A 

5330 

4150 

6160 

6010 

6460 

6480 

7770 

6640 

7950 

6720 

8220 

8230 

8450 

ii)  Aluminum,  Sodium,  and  Potassixim 

The  most  intense  lines  for  each  element  are  given, 
together  with  the  transition  involved,  the  transition  probability,  and  the 
energy  level  and  degeneracy  of  the  upper  state. 


Aluminum 


Transition 

X(i) 

A(sec 

g 

E(cm”^) 

3d  2d5/2^  3p 

^P3/2 

3092. 71 

6.95x1 0^ 

6 

32436. 79 

4s  ^Sl/Z  3p 

^P3/2 

3961. 53 

2. 21  X  10^ 

2 

25347. 69 

Sodium 

4p  ^P3/2  —3s 

^Sl/2 

3302. 32 

3. 78  X  10^ 

4 

30272. 51 

3p  2p3/2  ^3s 

^Sl/2 

5889. 95 

9. 04  X  10^ 

4 

16973. 38 

3d  ^D5/2  — 3p 

^P3/2 

8194. 81 

1. 39  X  10^ 

6 

29172. 86 
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Potassium 


6p  ^P3/2-^  4s 

3446. 72 

7.  32  X 

10 

4 

29007. 70 

5p  ^P3/2—  4s  ^81/2 

4044. 14 

3.  01  X 

10^ 

4 

24720. 20 

4p  2p3/2—  4s  2si/2 

7664. 91 

5.  57  X 

lo”^ 

4 

13042. 89 

iii)  Cesium 


The  first  table  is  similar  in  all  respects  to  the  pre- 
ceeding  one,  i.  e^  ,  transition,  wavelength,  probability,  degeneracy,  and 
energy  are  given.  The  second  table  presents  the  Stark  broadening  data  of 
Griem.  The  data  is  in  the  form  of  half-widths  of  the  specified  line 

(in  9i)  for  an  electron  density  of  10^^  cm“^  at  two  temperatures. 


Transition 

\  O 

X(A) 

A(sec 

g 

E(cm 

9p 

2p3/2—  6s 

^81/2 

3611. 52 

3.  49 

X 

10^ 

4 

27681. 96 

9p 

^Pl/2-»6s 

^81/2 

3617.41 

2.  35 

X 

10“^ 

2 

27637. 29 

8p 

^P3/2-*-  6s 

^81/2 

3876. 39 

1. 54 

X 

10^ 

4 

25791. 78 

8p 

^Pl/2—  6s 

^81/2 

3888. 65 

5.  52 

X 

10^ 

2 

25709. 14 

7p 

^P3/2-«-6s 

^81/2 

4555.  35 

4.  82 

X 

10^ 

4 

21946. 66 

7p 

^Pl/2  -►6s 

^81/2 

4593. 18 

2.  06 

X 

10^ 

2 

21765. 65 

8d 

^D3/2  — 6p 

^Pl/2 

6012. 14 

3.47 

X 

10^ 

4 

27811. 25 

8d 

^D5/2  -►  6p 

^P3/2 

6214. 81 

1. 54 

X 

10^ 

6 

27822. 94 

8d 

^D3/2  -►6p 

^P3/2 

6219. 33 

1. 71 

X 

5 

10 

4 

27811. 25 

7d 

^D3/2  -►  6p 

^Pl/2 

6725. 12 

1.42 

X 

10^ 

4 

26047. 86 

7d 

^D5/2  -►6p 

^P3/2 

6975. 21 

2.  63 

X 

10^ 

6 

26068. 83 

7d 

^D3/2  —  6p 

^P3/2 

6985.43 

2.  92 

X 

10^ 

4 

26047. 86 

6p 

^P3/2  -►  6s 

^Sl/2 

8521. 10 

4.  85 

X 

10^ 

4 

11732. 35 

6d 

^D3/2  ■—  6p 

^Pl/2 

8761. 38 

2.  80 

X 

10^ 

4 

22588. 89 

6p 

2pi/2  —6s 

2s1/2 

8943. 50 

2.  84 

X 

10^ 

2 

11178. 24 

6d 

^D5/2  -►  6p 

^P3/2 

9172. 24 

7.  00 

X 

10^ 

6 

22631.83 
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Transition 


o 

Half-Width  (A) 


2500  °K 

5000 

7p  ^P3/2 

6s^Sl/2 

0.  205 

0.  274 

7p  ^Pl/2 

—  6s  ^Sl/2 

0.  168 

0.  208 

8p  ^P3/2 

6s  ^Sl/2 

0,  613 

0.  718 

8p  ^Pl/2 

2 

6s  Sl/2 

0.  491 

0.  593 

9p^  P3/2 

6s  ^Sl/2 

1.  12 

1.  34 

9p  ^Pl/2 

—  6s  ^Sl/2 

1.21 

1.44 

7d  ^D3/2 

—  6p  ^Pl/2 

1. 30 

1. 63 

7d  ^D3/2 

—  6p  ^P3/2 

1. 40 

1. 76 

7d  ^D5/2 

6p  2p3/2 

1. 32 

1. 69 

8d  ^D3/2 

-►  6p  ^Pl/2 

3,  02 

3.  66 

8d  ^D3/2 

—  6p  ^P3/2 

3.  23 

3.  92 

8d  2d5/2 

—  6p  ^P3/2 

3.  05 

3.  77 

9p  ^P3/2 

—  5d  ^D3/2 

6.  73 

7.  70 
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iv)  Hydrogen 


This  section  presents  the  required  spectroscopic 
data,  including  Stark  Coefficients,  for  the  first  three  lines  of  the  Balmer 
Series  of  the  hydrogen  atom,  i.  e.  ,  Hq,  ,  ,  and  H-y. 

The  first-order  Stark  Coefficients  are  taken  from 
Condon  and  Shortley  and  are  given  in  Angstrom  units  per  kilo- 

volt/6m.  It  should  be  noted  here,  that  the  data  of  Condon  and  Shortley  do 
not  apply  to  the  usual  eigenfunctions,  characterized  by  the  quantum 
numbers  n,  1,  m,  in  a  spherical  coordinate  system.  Rather,  a  parabolic 
coordinate  system  is  employed,  characterized  by  the  quantum  numbers 
ki,  k2,  m. 


Here  n  =  kj^  +  k^  +  |mj  +  1  (J-1) 

and  n  and  m  have  the  same  significance  as  in  the  spherical  case.  Ac-*- 
cording  to  Schiff  (J“5)  ^  the  eigenfunctions  in  one  system  are  linear  com¬ 
binations  of  the  eigenfunctions  in  the  other  system. 

Thus,  since  the  Stark  data  is  given  for  transitions 
between  states  having  eigenfunctions  based  on  the  parabolic  coordinate 
system,  a  consistent  set  of  transition  probabilities  must  be  determined. 
The  transition  probability  is  given  by 


A 


nm 


64  € 


3hX- 


j^nm  j^mn 


2.  025  X  10^® 


(J~Z) 


where  S  is  the  line  strength  of  the  transition  expressed  in  units  of 
€  ^  a^^and  ^is  given  in  Angstrom  units.  Thus 
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s 


^2  ^nm 


Rmn 


+J’^n  'pm'  y  '^’’^'Pm  'Pn  V 

tj'/'n  'PJz  dT  ^\pm'P^'  ^  <3'^ 


Here  the  \//n  and  are  the  eigenfunctions  of  the  two  states  involved  in 

the  transition,  and  the  asterisk  indicates  the  complex  conjugate.  The  ex¬ 
pressions  for  X,  y,  z,  dTand  are  found  in  Condon  and  Shortley. 


o 

The  first  table  below  gives  the  wave  length  in  A,  and  upper  state 
energy  in  cm“  ^ ,  for  the  three  lines  of  the  Balmer  Series  considered. 

The  second  table  gives  the  transition  probability,  in  reciprocal  seconds, 
and  the  first-order  Stark  coefficient,  in  K  per  kilovolt/ cm,  for  each  specific 
transition  characterized  by  k^,  k^,  m. 


Line 

Wave  length 

Upper  State  Energy 

6562. 85 

97492. 

4861. 33 

102824. 

Hy 

4340.  47 

105292. 

line  (n  =  3  to  n  =  2) 

Transition 

Probability (x  10"^  ) 

Stark  Coefficient 

110-010 

.  0858 

.  0555 

101-001 

.  162 

.  0833 

200-100 

.  235 

.  1108 

002-001 

.  648 

0. 

110-001 

.  0124 

0. 

101-100 

.  272 

.  0277 
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line  (n  =  4  to  n 


2) 


Transition 

210-010 

201-001 

300- 100 
111-010 
102-001 
210-001 
201-100 

H  line  (n  =  5  to  n 

Transition 

220-010 

211-001 

310-010 

301- 001 
400-100 
112-001 
220-001 
211-100 
202-001 
310-001 
301-100 


Probability{x  lO'®) 

.  0663 
.  0513 
. 00962 
. 00480 
.  103 
.  0203 
. 00786 

2) 

Probability  (x  10“^) 
. 00548 
, 00337 
.  0694 
,  0109 
.  0462 
.  0403 
. 00928 
. 000399 
.  0269 
.  193 
.  0299 


Stark  Coefficient 
.  0914 
.  1218 
.  1523 
.  0304 
.  0609 
.  0609 
.  0914 


Stark  Coefficient 
.  0243 
.  0607 
.  1458 
.  1882 
.  2185 
0. 

0. 

.  0364 
.  1217 
.  1217 
.  1581 


B)  Molecular  Specie 

There  will  be  at  least  three  band  spectra,  due  to  molecular 
constituents  of  the  rocket  motor  exhaust,  appearing  in  the  wavelength 
region  of  interest.  The  higher  vibrational  transitions  of  the  first  posi¬ 
tive  band  of  nitrogen,  as  well  as  the  second  positive  band  of  nitrogen, 
and  a  band  in  the  near  infrared  due  to  oxygen  should  be  encountered. 
The  data  below  give  the  transitions  involved,  the  locations  of  the  band 
heads,  and  the  molecular  spectroscopic  constants. 
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Nitrogen 


First  Positive 

Xoo  =  10470  A 


Second  Positive 
X  oo  =  3370  A 
c  B  ^TT 

-  -U  ^  'g 


State 

T(cm"^)  CO 

e(cm"^)  a)eXg(cm"^ 

J.  g(cm-^) 

89147. 

2035.  1 

17.  08 

1. 8259 

0. 0197 

B  ^r]- 
'  g 

59626. 

1734.  11 

14.47 

1. 6380 

0. 0184 

A  ^X 

50206.  0 

1460. 37 

13. 891 

1. 440 

,  0.013 

Oxygen 

X  oo  =  7620  A 
b  -  X 

g  ^-g 

State  T(cm  OJ  (cm"^)  COx  (cm~^)  B„(cm"^)  CL  (cm"^) 

^  0  0  G  0  ^ 

1  + 

b  Yg  13195.222  1432.687  13.9500  1.40041  0.01817 

X  ^Xg  0  1580.361  12.0730  1.44566  0.01579 
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